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We study the contribution of polyelectrolytes 
in solution to the bending moduh of charged 
membranes. Using the Helfrich free energy, 
and within the mean-field theory, we calcu- 
late the dependence of the bending moduli 
on the electrostatics and short-range interac- 
tions between the membrane and the polyelec- 
trolyte chains. The most significant effect is 
seen for strong short-range interactions and low 
amounts of added salt where a substantial in- 
crease in the bending moduli of order 1 ksT 
is obtained. From short-range repulsive mem- 
branes, the polyelectrolyte contribution to the 
bending moduli is small, of order 0.1 ksT up to 
at most 1 ksT. For weak short-range attraction, 
the increase in membrane rigidity is smaller and 
of less significance. It may even become neg- 
ative for large enough amounts of added salt. 
Our numerical results are obtained by solving 
the adsorption problem in spherical and cylin- 
drical geometries. In some cases the bending 
moduli are shown to follow simple scaling laws. 



1 Introduction 

The study of interactions between charged and 
fiexible membranes and polyelectrolytes (PEs) 
in solution has generated a lot of interest in re- 
cent years, partly motivated by the importance 
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of such interactions in biology. Understand- 
ing the interaction between charged macro- 
molecules such as DNA, RNA and various pro- 
teins and biological cell membranes (modeled 
as charged and fiexible interfaces) sheds light 
on many important cellular processes. Besides 
its biological significance, the adsorption of PEs 
onto charged membranes raises interesting ques- 
tions about the interplay between short-range 
and electrostatic (long range) interactions in 
these multi-component charged systems. Re- 
cent works include numerical and analytical cal- 
culations [IlISllsllllElEllTllHlElIToliniand 

scaling arguments PQIH El [121 [13] • 

In our previous work [14J, three regimes 
have been found for polyelectrolyte adsorption, 
(i) When the short-range interaction between 
the membrane and the PE chains is repulsive, 
the surface charge is low and the ionic strength 
of the solution is high, the polymers deplete 
from the charged membrane, (ii) For higher sur- 
face charges, or lower ionic strength, the poly- 
electrolytes adsorb on the charged membrane 
and screen the surface charges, (iii) When the 
short-range interactions between the membrane 
and the PE chains are attractive, the PE chains 
adsorb on the membrane, and the adsorbed 
layer carries a higher charge than that of the 
bare membrane. In this situation the polyelec- 
trolytes over compensate the bare surface — 
a phenomenon of practical importance in the 
build-up of multilayers of alternating cation and 
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anion polyelectrolytes 

The flexibility of fluid surfaces and membranes 
has been studied in various cases, and it depends 
on the lipid composition, tail length and molec- 
ular tilt. In the case of charged and flexible 
membranes immersed in a pure ionic solution 
(no macromolecules) [161 [IH [HI [13 [201 [2ll [22l 
[231 [21 [251 [261 [27], the electrostatic contribu- 
tion was found to increase the membrane rigid- 
ity, and the addition of salt to decrease it. In a 
different set of studies, the adsorption of neutral 
polymers on membranes has been investigated 
theoretically [281 [291 [30l [3ll [32l [331 [M] and in 
experiments [351 ISSl [371 [38] . It was found that 
the addition of polymers reduces the membrane 
rigidity. In all those cases, the contribution of 
the solution to the bending rigidity was found to 
be of order 0.1 — 1 fc^T, which is low in compar- 
ison to the intrinsic monolayer bending rigidity 
of approximately 10 fc^T. We note that in spe- 
cial cases, by adding a co-surfactant (alcohol), 
the membrane bending rigidity can be brought 
down to values of roughly 1 ksT |39j . 

In the present work, we study the combined 
system of charged chains interacting with op- 
positely charged and flexible membranes. Our 
study is similar in spirit to that done for DNA- 
lipid systems [^ [JT] , where the DNA was mod- 
eled as rigid and charged rod. The main dif- 
ference is that our charged macromolecules are 
flexible. The above mentioned three regimes 
dictate a different contribution to the membrane 
rigidity and stability. For short-range repulsive 
membranes, the membrane rigidifles due to its 
charges while its rigidity decreases with the in- 
crease of the polyelectrolyte charge. For weak 
short-range interactions, the contribution to the 
membrane rigidity decreases and may become 
slightly negative, while for strong short-range 
attraction the membrane becomes rigid again. 
In most cases, the magnitude of the contribu- 
tion to the bending rigidity is of the same or- 
der of magnitude as that of neutral polymers or 
pure salt solutions, namely 0.1 ksT up to 1 /c_bT. 
However, we show in this paper that for strong 
enough short-range attraction a more significant 
contribution, of order ksT or higher can be ob- 
tained. 



2 Mean Field Equations 

2.1 Free-Energy Formulation 

Consider a bulk aqueous solution containing 
polyelectrolyte (PE) chains, along with their 
counter ions and added salt. A curved and 
charged membrane is placed at the origin |r| = 
0. In order to extract the membrane elastic 
moduli, we take the membrane shape to be ei- 
ther spherical or cylindrical with radius R . The 
free energy for the combined system has been 
formulated before [T], and can be written as a 
sum of four contributions Ftot = Rpoi + -^ions + 

-^el + -^int- 

The polymer free energy, Fpoi, is: 
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2 is the monomer 



where a is the monomer size, 
bulk concentration, (j)"^ is the local monomer 
concentration. The free energy is given in terms 
of ksT, and the integration is carried out over 
the entire volume outside the sphere/cylinder, 
|r| > R. The first term in Eq. ([T|) accounts for 
the chain elasticity. The second and third terms 
account for the excluded volume interactions in 
the solution and in the bulk, respectively: 

E (0) = (1 - a302) log (1 - a302^ _ 

- a'<P') <P' (2) 
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where v is the excluded volume coefficient. For 
0^ <^ the above expression can be expanded 
to the well known |f (0^ — 0^) , which is com- 
monly used for the excluded volume interaction. 
The small ion contribution to the entropy is: 
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where c± and are the local and bulk concen- 
trations of the ± small ions, respectively. The 
third contribution to the free energy, the elec- 
trostatic free-energy, is: 



dr (c+ - c_ + C - 
-i- / dr|VC|'+ dAaC(5) 

SttIb JV J\r\=R 

where ( = eip/ksT is the renormalized elec- 
trostatic potential, / is the fraction of charged 
monomers and Ib = e^/eksT is the Bjerrum 
length. The first integral accounts for the in- 
teractions between the positive ions, negative 
small ions and the monomer charges with the 
electrostatic potential. The second integral ac- 
counts for the self-energy of the electrostatic 
field and the third to the interaction between 
the surface charge and the electrostatic poten- 
tial. Note that the third integral is taken only 
over the charged body surface, |r| = R. 
The last part of the free energy is the short- 
range interaction of the PE chains with the 
membrane: 



int 



dA, 
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6d J\r\=R 

where (j)s = (j){\r\ = R) is the monomer concen- 
tration on the surface. This term stands for a 
general short-range interaction, where the in- 
teraction strength is determined by the phe- 
nomenological constant d~^, which has dimen- 
sions of inverse length. 

Minimization of the total free energy Ftot = 
-^poi + -Pions + -^ei + -^int yields the following mean- 
field equations, expressed in terms of a dimen- 
sionless variable rj = 
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where = (SvrZ^Csait) is the Debye-Hiickel 
length scale for the screening of the electro- 
static potential in the presence of added salt 



and = {47ilB(plf) ^^"^ the corresponding 
length for the potential decay due to counteri- 
ons. Note that the actual decay of the electro- 
static potential is determined by a combination 
of salt, counterions, and polymer screening ef- 
fects. The excluded volume interaction is taken 
into account using the function: 



A{V) = log(l - 0^a^) - log(l - 0,Vr/^) + 
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(9) 



which represents the full excluded volume inter- 
action. 

The solution of Eqs. ([7]) and ([8]) requires four 
boundary conditions. Two of the boundary 
conditions are taken far from the membrane, 
where the monomer concentration and elec- 
trostatic potential retrieve their bulk values 
?7(|r| oo) 1 and C(kl ~^ oo) 0. The 
other two boundary conditions account for the 
interaction with the charged membrane, and 
can be obtained by integrating Eqs (I7j) and (E]) 
from |r| = i? to a small distance from the mem- 
brane, yielding: 



A- VC||r|=R 



—AtvIbct 



(10) 
(11) 



Equation ffTOl) is the usual electrostatic bound- 
ary condition for a given surface charge density, 
while Eq. ( fTTl) is the Cahn - de Gennes boundary 
condition [42j, which is often used for calculat- 
ing polymer profiles [30l [311 [29] . 
For large i?, the total free-energy can be ex- 
panded around its fiat surface value in the fol- 
lowing way [B]: 
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where Ftot is the total free energy (in units 
of kBT), /o is the free energy per unit area 
of a solution in contact with a planar surface 
{R —>■ oo) that has the same system parameters. 



3 



and Ci , C2 are the radii of curvature for the mem- 
brane. For spherical surfaces the radii of cur- 
vature are ci = C2 = 1/-R, while for a cylindrical 
surface Ci = C2 = 0. The parameters 5k and 
5hg are the contributions of the PE (and salt) 
solution to the mean and Gaussian curvature 
moduli of the membrane, respectively, in units 
of UbT. Namely, n = k^ + 6k and kg = '^g + '^'^g' 
include the intrinsic values as well as the contri- 
butions coming from the solution. Throughout 
this paper we will consider only the changes in 
the elastic moduli k, kg with respect to their 
bare values. An increase in the mean curva- 
ture modulus 5k increases the membrane rigid- 
ity, while an increase in the Gaussian modulus 
5kg makes saddle points on the membrane more 
favorable. 

The parameter Rq is the radius of spontaneous 
curvature for the membrane, which is depen- 
dent on the exact chemical composition of the 
membrane. For a positive Rq, the membrane 
bends towards the external solution, while for 
negative Rq it bends away from the solution. 
In the case of a bilayer membrane, where the 
same solution is in contact with both leaflets of 
the membrane, the spontaneous curvature van- 
ishes due to symmetry. In this paper, we do not 
calculate the spontaneous curvature but rather 
focus on the changes to the bending moduli 5k 
and 5kg- 

2.2 Numerical Procedure 

Equations (17|)- (fTT|) are solved numerically for 
the cases of a charged sphere and a charged 
cylinder. The numerical procedure follows the 
relaxation scheme [13] , as was described in pre- 
vious publications [HE]. For each solution, we 
calculate the total free energy per unit area 
for the cases of a spherical membrane fs and 
a cylindrical membrane fc- The contributions 
of the polyelectrolyte solution to the mean and 
Gaussian curvatures are then calculated by ex- 
panding /c, fs to second order in 1/R: 

fs = /o + ^ + |l (13) 
fc = /o + 42 + ^ (14) 



The contributions of the solution to the curva- 
ture moduli are given by: 

5k = 2Bc (15) 
5kg = Bs-^Bg. (16) 

A surface is stable under long wave bending 
fluctuations only when k > 0, and against spon- 
taneous vesiculation (topological change) when 
2k + Kg > [21]. In the following we show that 
the contribution of the PE solution to the stabil- 
ity of a charged surface depends on the amount 
of added salt, and has a non-monotonic depen- 
dence on the short-range interactions between 
the membrane and the polyelectrolyte. 

3 Results 

We find large contributions, of order 1 ksT, 
to the surface bending moduli for the case 
of strong short-range attractive surfaces. For 
weaker short-range interactions and for repul- 
sive surfaces, the contribution is smaller, of or- 
der 0.1 — IksT. We discuss first the strong 
repulsive and strong attractive surface limits, 
and then turn to the intermediate case, where 
the contribution to the bending rigidity is less 
significant. We present analogies and scaling 
calculations to explain the different regimes. 

3.1 Strong Short-Range Attrac- 
tive Membranes 

In previous publications [5l [6], the adsorbed 
amount of PEs as well as the layer width were 
studied in detail for short-range attractive sur- 
faces. The adsorbed PEs charge was shown 
to exceed the bare surface charge significantly 
for strong short-range interactions, and to ex- 
ceed it mildly for weak short-range interactions. 
The width of the adsorbed layer, on the other 
hand, depends on the shorter of the two ad- 
sorption length scales: (i) d, the length scale for 
short-range attraction from Eq. flTT]) . and (ii) 

^ = a/ \fv^, the Edwards correlation length for 
neutral polymer adsorption. In our model, we 
use an almost theta solvent < vcj)^ ^ 1, so we 
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assume that the shorter length scale is always 
c? -C ^. The increase in short-range attractive 
interactions, in this case, decreases the adsorbed 
layer width [4j, [5] . The combined charge of the 
PE-membrane complex is, therefore, opposite 
to the initial surface charge, and its magnitude 
may be much higher than the initial membrane 
charge. This charge is distributed within a layer 
of width d close to the surface. 
Our focus in this section is on the low salt case. 
In Fig. 1 we show the dependence of k and 
on the short-range interaction parameter 
for low salt conditions. For strong short-range 
attraction, namely d < 2A, the calculated 5n 
is positive and increases strongly, reaching val- 
ues of several ksT. The 5kg is negative, and 
shows a stronger increase with d~^ than the cor- 
responding 5n. The contribution to the vesicu- 
lation stability 25k+5kg in this case can be seen 
to be negative for high c?"^, leading to destabi- 
lization of the surface. For lower (but positive) 
d~^, the contribution to 25 n + 5kg is positive, 
and thus enhances the membrane stability. 
These results can be explained by the following 
argument. When the adsorbed layer width is 
smaller than the electrostatic screening length, 
the membrane-PE complex can be viewed as a 
renormalized charged surface, containing both 
the bare surface charges and the adsorbed PE 
charges, in a weak ionic solution. The renor- 
malized surface interacts with the ionic so- 
lution in the same manner as a bare mem- 
brane Uni m [Tg [m [201 [211 [221 [231 [2l|. The 
strong increase in the surface charge and the 
lack of small ion screening (low salt) cause the 
surface fluctuations to be strongly unfavorable, 
making 5n positive. This allows substantial in- 
crease in the magnitude of (5k, amount- 
ing to several fc^T, which is a significant con- 
tribution to the membrane curvature moduli, 
as can be seen in Fig. 1. We note that as 
the renormalized surface charge increases, the 
value oi 5k should approach the low salt limit 
for charged surfaces in pure ionic solutions (no 
PE) 5k XD/i^nlB) ~ TfceT [21] (See Fig. 1). 
However, this limit is still higher than our nu- 
merical results for the PE-membrane complex. 
The crossover between positive and negative 
values of 25k + 5kg can also be explained 



by analogy to ionic solutions. There are two 
regimes for 25k+5kg in pure ionic solutions [23], 
depending on the surface charge. For weakly 
charged surfaces Csait ^ (or \d <^ 

Agc where Aqc = (2vr/B is the Gouy- 

Chapman length), 5 well as 25k + 5kg 
are positive [181 120], while for highly charged 
surfaces Cgait -C lB\cr\ we get 5k > and 
25k + 5kg < [211 [22] ■ In our case of poly- 
electrolyte solutions, for low enough attractive 
short-range interactions the renormalized sur- 
face charge is low, and 25k + 5kg > 0. For 
stronger short-range interactions, the renormal- 
ized surface charge increases and 25k + 5kg < 0. 
The results in Fig. 1 were presented in the low 
added salt case, d <^ Xd- When more salt is 
added into the solution the electrostatic interac- 
tions between the polymers become weaker. In 
case of a high amount of added salt, d^ Xd, the 
charges of the adsorbed polymers are screened 
over smaller length scales than the adsorbed 
layer width, and the analogy to a renormalized 
charged surface breaks down. We further dis- 
cuss this case in Sec. 13.31 

3.2 Short-Range Repulsive 
Membranes 

Without the PE adsorption, the charges on the 
membrane increase the membrane rigidity, k, 
due to the long-range repulsion between them. 
For strongly repulsive surfaces, the adsorption 
of PE chains to the membrane screens its surface 
charges and causes k to decrease. Note that the 
PE charges do not fully compensate the mem- 
brane bare charges like they do in the strongly 
attractive regime, leading to the difference in 
the corresponding system behaviors. 
There are two regimes for the case of short-range 
repulsive membranes: (i) the low-salt regime 
where the surface charges are mainly balanced 
by the adsorbed monomer charges, and (ii) the 
high-salt regime where the salt ions screen the 
surface charges and the monomers deplete. The 
bending moduli in the latter regime are simi- 
lar to those of a strong ionic solution with no 
added PE, as were discussed in detail in other 
papers [161 [IH [HI [201 [21] . 
In Ref. [5 we addressed the PE adsorption 
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close to a flat charged surface and showed that 
the screening length of the electrostatic poten- 
tial depends strongly on the amount of added 
salt. For low salt conditions, the surface charges 
are screened mainly by the adsorbed monomer 
charges. In this case the free-energy per 
unit area scales like /ads ~ {af^"^ f~^^^a'^^^f^^ 
and the screening length scales like D ~ 
a?/^ /{fls For high amounts of added 

salt, the screening is mainly done by small ions 
and the PE chains deplete. The transition from 
depletion to adsorption regimes occurs when the 
screening length due to monomer adsorption be- 
comes similar to the one for small ion adsorp- 
tion, i.e., when A/) ~ -D. 

The fiat surface results can now be easily ex- 
tended to curved surfaces. In the low-salt 
regime, monomers are adsorbed to the mem- 
brane, and the length scale for the adsorp- 
tion is D. We expect the free energy of the 
curved membrane (per unit area) to scale like 
/tot = fadshiD/R). Expanding /tot to second 
order in and comparing to Eq. (1121) shows 
that both curvature moduli scale like: 

I I 2 

5n, 5kg ~ UsD' ~ (17) 

which is of order 0.1 — 1 ksT for physiological 
range of system parameters. 
In Fig. 2 we present the values of 6n and 6kg 
as a function of the surface charge \a\, in the 
case of strongly repulsive membranes = 
— 20A~^ and monomer size of a = 10 A. For 
low amounts of added salt we find a scaling re- 
lation Sk, 5 Kg ~ |ct|^ with j3 ~ 1.2. Note that 
the numerically calculated exponent /? ~ 1.2 is 
slightly larger than (3=1 derived in Eq. (ITTl) . 
This discrepancy seems to occurs because the 
high monomer size used here makes the full ex- 
cluded volume interaction substantial. For com- 
parison, we plot the corresponding 5k, 5kg for 
the case of a pure ionic solution, namely when 
no polyelectrolytes are added to the solution. 
In this case the magnitudes of both 5k, 5kg are 
very high. This shows that the addition of poly- 
electrolytes into low salt solutions can cause a 
strong reduction of the bending moduli, in the 
order of several ksT as compared with the pure 
salt solution. This reduction can be explained 



by the fact that the polyelectrolytes replace the 
salt ions in screening the surface charges, thus 
allowing greater membrane flexibility. 
In the high salt regime, the membrane interacts 
only with the small ions, and the polymer chains 
are depleted. The free energy in this case is 
similar to that of a weakly charged surface in 
an ionic solution [23]. Namely, the screening 
length is A^i and the free energy per unit area 
for the case of a fiat surface is /dcp ~ \a\^lB)^D- 
The curvature moduli are the same as in a pure 
ionic solution [HI [HI [Tg [201 [23] , where: 

5k, 5kg ~ /depA|, ~ |<t|^ IbX^d- (18) 

Note that the depletion condition Xa < D im- 
plies that the bending moduli for the high salt 
(depletion) case are always lower than for low 
salt (adsorption) case, despite the depletion of 
the polymers. Both 5k and 5kg in this case are 
of order 0.01 — 0.1 ksT for physiological range 
of system parameters, and scale like |(t|^. 
For both low and high salt regimes, the contri- 
bution to the Gaussian bending modulus 5kg is 
negative, resulting from the electrostatic repul- 
sion between the membrane constituents |19j . 
The contribution of the PE solution to the vesic- 
ulation stability 25k + 5kg is always positive, 
even for very low salt concentrations, in contrast 
to the pure ionic solution results [2T1 [221 [25] . 
The low-salt regime of the ionic solutions, in 
which the contribution to 25k + 5kg is nega- 
tive, is replaced here by the adsorption regime, 
in which this contribution is still positive. 

3.3 Weak Short-Range Interact- 
ing Membranes 

In Fig. 3 we present an enlargement of Fig. 1 for 
surfaces having only a weak short-range inter- 
action with the PE in solution, ~ 0. As can 
be seen, the magnitudes of both 5k and 5kg 
decrease substantially for low and may 

become negative for high values of added salt. 
The decrease in the magnitude of the bending 
moduli can be attributed to the strong screen- 
ing of the surface charges by the adsorbed PEs, 
which makes 5k smaller. We note that the con- 
tribution to 5k and 5kg is negligible in compar- 
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ison to the membrane intrinsic bending moduh, 
and probably cannot be observed experimen- 
tally. This regime is presented only in order 
to complete the dependence picture. 

For low and positive and high amounts of 
added salt, we find a negative contribution to 
the mean curvature modulus due to the poly- 
mer adsorption 6k < 0, as well as in 26k + 6kg- 
This surprising result can be explained by anal- 
ogy to the neutral polymer case. In past pub- 
lications [221 EDI EI] it was shown that neu- 
tral polymer solutions have negative 6k, posi- 
tive 6kg and negative 26k + 6kg. We find here 
similar results. 

The analogy to neutral polymer adsorption is 
important and can be understood in the follow- 
ing way. The high amount of added salt screens 
both the surface and the PE charges, so their 
effective interaction becomes short-ranged lead- 
ing the way to an almost neutral polymer behav- 
ior. The increase in 6kg derives from the effec- 
tive attraction between membrane constituents, 
which results from their short-range attraction 
to the PE chains. The analogy to neutral poly- 
mers requires the screening length for the elec- 
trostatic interactions to be much smaller than 
the layer width Xd < d. This is satisfied in 
high salt and low short-range attraction con- 
ditions. For stronger short-range interactions, 
the layer width decreases, and the electrostatic 
interactions between monomers becomes impor- 
tant. In this case, the analogy to neutral poly- 
mers breaks down, and 6k, 26k + 6kg start to 
increase back. For higher d~^, both 6k and 
26k + 6kg become positive again, marking the 
crossover to the strong attraction regime de- 
scribed in Sec. 13.11 above. This can be seen in 
Fig. 3, where for high amounts of added salt, 
the increase in 6k with d~^ indeed begins when 
A/5 ^ (i, as expected from the above analysis. It 
is also important to note, that the magnitude 
of both 6k and 6kg are very small, of order 
0.01 ksT, and are negligible in comparison to 
the intrinsic bending moduli of a membrane, of 
order 10 ksT. The decrease shown here cannot 
account for a decrease in membrane rigidity that 
was recently found for DNA adsorption 



4 Conclusions 



In this paper we show that the interaction be- 
tween charged and fiexible polymers (PE) and 
oppositely charged and fiexible membrane de- 
pends both on their electrostatic and short- 
range interactions. A non-monotonic depen- 
dence of the curvature moduli is obtained as 
function of the short-range interaction between 
the membrane and the PE chains. We find a 
significant contribution to the bending moduli, 
of order of several ksT, in the case of strong 
short-range attraction between the PE chains 
and the surface. For weak attractive interac- 
tions, the contribution of the PE solution to the 
membrane curvature moduli is small (in units of 
ksT), and for repulsive interactions it increases 
back, and may reach values of 0.1 — 1 A;_bT. 

Our work deals only with uniformly charged 
membranes. In biological membranes, however, 
the membrane is composed of a mixture of neu- 
tral and charged lipid. The lipid molecules 
can rearrange and can cluster around oppositely 
charged PE [17] . This in turn can have a strong 
effect on the overall membrane rigidity. Future 
works may offer extensions of the present one by 
calculating the contributions to the spontaneous 
radius of curvature Rq, especially in the case 
of asymmetrical solutions, or when the mem- 
branes are composed of asymmetric inner and 
outer leaflets. Other potential directions may 
include changes in the effective lipid headgroup 
size and water activity due to the presence of 
polyelectrolytes. 
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Figure 1: The d ^ dependence of (a) 5k and (b) Skq is presented for low added salt concentration 
(csait = 0.1 mM). Other parameters are \a\ = 0.001 A"^ a = 5 A, / = 0.5, v = 50A'l 0g = 
lO-^A-3, T = 300K and e = 80. For large d-\ we see a significant increase in the magnitude of 
both 6k and Skq, amounting to several ksT, which is very significant in comparison to normal 
membrane curvature moduli. 
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Figure 2: The dependence of 5k (a) and Skq (b) on the surface charge is presented for the case of 
repulsive membranes. The triangular symbols are the numerically calculated 5k for Csait = 0.1 mM, 
a = lOA, / = 0.5, V = 50 A^, 02 = 10-8 A-3, T = 300K and e = 80. The solid line scales as 
\a\^ with P ~ 1.2. As can be seen, for low salt concentrations the exponent of the curvature 
modulus, P ~ 1.2, is close to the predicted (3 = 1 derived from Eq. flTTI) . The dashed-dotted hne 
is the numerically calculated 5k, 5kg for the case of an ionic solution with no polymers, with 
parameters Csait = 0.1 mM, T = 300 K and e = 80. The addition of polyelectrolytes can, in this 
case, reduce the curvature moduli significantly. 
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Figure 3: The dependence of 5k and Skq is presented for two salt concentrations: (a) 
Csait = 0.1 M and (b) Csait = 0.1 mM. Other parameters are \a\ = 0.001 A~^, a — 5 A, / = 0.5, 
V = 50 A'"^, (t)l = 10^^ A^^, T = 300 K and e = 80. In both plots, the solid line corresponds to 6k 
and the dashed one to Skg- Three regimes for both curvature moduli are seen. For < we 
obtain Sk > and 6kg < 0, for ~ the sign of both is inverted, and for » both moduh 
return to their original sign. The magnitude of both moduli is very small for this parameter 
range. 
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